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Abstract 

A generalized optimal velocity model is analyzed, where the optimal velocity 
function depends not only on the headway of each car but also the headway 
of the immediately preceding one. The stability condition of the model is 
derived by considering a small perturbation around the homogeneous flow 
solution. The effect of the generalized optimal velocity function is also con- 
firmed with numerical simulation, by examining the hysteresis loop in the 
headway-velocity phase space, and the relation between flow and density of 
cars. In the model with a specific parameter choice, it is found that an in- 
termediate state appears for the movement of cars, where the car keeps a 
certain velocity though the headway is short or long. This phenomenon is 
different from the ordinary stop-and-go state. 
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1 Introduction 



Traffic flow problem has been extensively studied from physical point of 
view. Fluid-dynamical model [|l|, cellular automaton model 0, and the car- 
following model0 have been proposed and analyzed in detail to understand 
the mechanism of the traffic congestion on a freeway. 

Toward a realistic model which explains the traffic flow dynamics, the op- 
timal velocity (OV) model proposed by Bando, Hasebe, Nakayama, Shibata, 
and Sugiyama[|, Q has attracted considerable interest. Based on the second- 
order differential equations, the model reveals the density pattern formation 
of the congested flow of traffic without introducing a time lag caused by the 
driver's response. 

Although the OV model is shown to have a universal structure in spatio- 
temporal patterns in the congestion, most of the analyses of the model have 
been done in the case where the optimal velocity function depends only on 
the headway of each car. One of the approaches to generalize the OV model 
is that the backward reference function is introduced^. Another approach 
to extend the OV model is to take into account the next-nearest-neighbor 
interaction 0, where the optimal velocity function depends not only on the 
headway of each car but also on the headway of the immediately preceding 
one. The generalized optimal velocity function is determined by taking into 
account the driver's skill, experience, and psychological effect, so that it is 
expected to describe more realistic traffic flow. 

The purpose of this paper is to analyze the generalized optimal velocity 
model proposed by Nagatani[0, |]. This paper is organized as follows. In 
Section 2 the generalized optimal velocity model is reviewed and its stability 
condition around the homogeneous flow solution is derived without long- 
wavelength approximation. In Section 3 numerical simulations are carried 
out, in particular, the hysteresis loop in the phase space and the flow-density 
relation are examined. In Section 4 we reconsider our model to compare the 
result with the one in the original model. In the model with a speciflc param- 
eter choice, it is found that an intermediate state appears for the movement 
of cars, which is different from the ordinary stop-and-go state. 

2 Generalized Optimal Velocity Model 

We flrst consider a dynamical model of the traffic flow given by 

Xn{t) = a(V(Ax„(t), AXn+l{t)) - Xn{t)), (1) 
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where Xnit) is the position of the n-th car at time t, Ax„(t) = Xn+iit) — x„(t) 
represents the headway of the car, and a is the sensitivity. Thus Axn+i(t) 
is the headway of the immediately preceding car. Here n = 1,2, ■ ■ ■ , N is 
each car number with N being the total number of vehicles. The driver's 
sensitivity a is assumed to be independent of n. Throughout this paper, we 
will consider the periodic boundary condition with respect to the coordinate 
Xn with period L. 

At first, let us look for an appropriate form of the optimal velocity func- 
tion to be suitable for our purpose. The driver sometimes pays attention to 
not only the headway but also the headway of the immediately preceding 
one. If the headway of the preceding car is short, the driver assumes that 
the forward driver decelerates, thus the driver decrease the optimal velocity 
even though the headway of his car is long enough. On the other hand, if the 
headway of the preceding car is long, the driver assumes that the forward 
driver accelerates, thus the driver increase the optimal velocity even though 
the headway of his car is short. 

Let us look at Figure 1, which describes the original optimal velocity 
function in Axn, Axn+i, and V space, where the numerical values of axes 
are not important here. The appropriate form of the function which satisfies 
the above requirement will be seen in Figure 2. 




Figure 1: The original optimal ve- Figure 2: The generalized optimal 
locity function V{Axn)- velocity function V{Axn, Ax„+i). 

As is already introduced in Ref.|^, |^, we can adopt the generalized opti- 
mal velocity function of the form 

V{AXn, AXn+l) = (1 - p)V{AXn) + pV{AXn+l), (2) 

where p is assumed to be independent of n and satisfies ^ p < ^ because 
the dominant part of the optimal velocity function should be Axn dependent 



3 



term. The advantage of the above form is to be able to find the effect of 
the additional term with varying p. In spite of introducing p and Aa;„+i 
dependent term, the above model provides exactly the same homogeneous 
flow solution as the original one without p dependence. When p = 0, the 
model reduces to the original one. According to the original OV function 
proposed by Bando et al.0], we will take hyperbolic tangent function of the 
form 

V{Axn) = tanh(Ax„ - 2) + tanh(2). (3) 

As we can see from eq.(H) with eq.(§), the generalized optimal velocity func- 
tion has the required properties. In general, p might depend on n and also 
depend on time t. The numerical simulation can be still performed under 
these circumstances. 

More generally, we can consider the form 

V{AXn, AXn+l) = U (Ax„) + W{AXn+l), (4) 

where U and W have the following properties: (i) monotonically increasing 
functions, (ii) they have upper bounds, and (iii) they satisfy \U\ > \W\ as a 
realistic model. 

Now, let us analyze the generalized OV model given by eq.(|I]) with eq.(^ 

Xn{t) = a {U{AXn{t)) + W{AXn+l{t)) - Xn{t)) ■ (5) 

As is easily seen from eq.(|^), they have a homogeneous flow solution 

x^^\t) = bn + ct, (6) 

where b = L/N and c = U{b) + W{b). We examine the stability against a 
small perturbation ?/„(t) around the homogeneous flow solution (j^). Substi- 
tuting 

Xnit)=xi^\t)+y^{t) (7) 
into eq.(^, the linearized equation to y„(t) is obtained as 

Vnit) = a {gAynit) + /iA?/„+i(t) - y^it)) , (8) 

where A?/„(t) = y^+iit) — ynit), and g and h is the derivative of U and W at 
b, respectively. The solutions to eq.(|D is given by the Fourier series as 

l/n,fc(t) = exp(iafcn + zt), (9) 

where = 2Tik/N with /c = 1, 2, ■ ■ ■ , and z satisfies 

z^ + az-a (^(e*"* - 1) + /^(e^*"^ - e^°^)) = 0. (10) 
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The stability condition is to find the condition Kez < for all modes a^. 
After a short calculation, it is equivalent to find the relation 

F{Y) = 32h'^Y^ + 16h{g-3h)Y'^ + 2 {{g - 3hf - 2ah) Y-a{g-h) < (11) 

holds for all modes a^, where we put Y = cos^ 

In order to proceed the analysis of stability condition, we will take the 
optimal velocity function adopted in eq.(||). Eq.([TI|) is rewritten by 

F{Y) = 32p'^fY^ + lQp{l-Ap)fY^+2{{l-4pyf-2apf)Y-a{l-2p)f < 0, 

(12) 

where / = V'{b). We can finally find that if the condition 

/<|(l + 2p) and p^i (13) 

is satisfied, eq.(0) holds for all modes ak- It should be emphasized that the 
stability condition (0) is derived without long-wavelength approximation. 
In the long- wavelength approximation, only the former condition in eq.(|l3D 
is derived from the stability analysis. The uniform solution is unstable if 
/ > f (1 + 2p) or p > |. Comparing the result with the original OV model, 
we can conclude that the model is stabilized in the region | ^ / < |(1 + 2p) 
by the effect of introducing the headway of the preceding car. 

3 Numerical Simulations 

To convince the analysis of stability condition for the generalized OV model, 
we will now solve 

Xn{t) = a ((1 - p)ViAx„) + pViAxn+i) - Xnit)) (14) 

numerically. In the simulation, a = 1 is taken throughout this paper. As 
the density of cars (p = N/L) varies and if / > f (1 + 2p) is satisfied, a 
homogeneous flow becomes unstable and makes a phase transition from free 
flow to congested one. Looking at the spatio-temporal pattern, the congested 
patters propagate backward. These characteristic features in the generalized 
model are the same as those appear in the original OV model. 

One of the typical features of the model is that the movement of the car 
becomes the stop-and-go states in the congested region and the congested 
pattern is very stable. It is well understood by examining the hysteresis 
loop in the headway-velocity phase space. In numerical simulations, we take 
N = 100 and L = 200 as an example, where / takes maximum value, because 



5 



the congested region is of our interest and the result depends only on the 
density but not on the number of cars and the circuit length. The initial 
condition we considered here is the homogeneous flow with small fluctuation, 
i.e. 

XniO) = bn + y^{0),Xn{0)=c, (15) 

where 2/„(0) is taken to be an uniform random distribution between —0.5 and 
0.5. 

2 
1.5 

> 

■B 1 

o 

0) 

> 

0.5 


0.5 1 1.5 2 2.5 3 3.5 4 
headway A 

Figure 3: The orbit of a car in the headway-velocity phase space after the 
organization of the congestion for p = and p = 0.2. The dashed curved line 
denotes the optimal velocity function. 




Figure ^ shows the orbit of a particular car in the head way- velocity phase 
space with the parameter p = 0, 0.2 as an example. After about 1,000 time 
when the generation of the congestion is finished, the shape of the hysteresis 
loop obtained in Figure 3 never changes. Furthermore, the shape does not 
depend on the initial random distributions. 

The numerical result shows that the effect of increasing p seems to be 
equivalent to increasing the sensitivity a. Of course, this is also guessed 
from the analysis of the stability condition in eg . (|I3|) . However, as will be 
seen in Section 4, we will find that the change of the value of p can not be 
compensated by rescaling the sensitivity a. 

Examining the bottom end point (Axc,fc) = (Aa;„,x„) and the top end 
point {Axi,Vf) = {Axn,Xn) in the hysteresis loop, we can obtain the back- 
ward velocity of the congestion given by 

ViAxc-VcAxf 
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Numerical simulation shows that the backward velocity of the congestion in- 
creases as p increases. We have also confirmed by simulating spatio-temporal 
patterns. The numerical results are listed in Table 1. 



p 


Axc 








Kack 


0.0 


0.32274 


0.03152 


3.67726 


1.89653 


0.14791 


0.1 


0.62051 


0.08319 


3.37945 


1.84485 


0.31302 


0.2 


0.91196 


0.16787 


3.08804 


1.76019 


0.49945 


0.3 


1.18567 


0.29206 


2.81434 


1.63600 


0.68632 


0.4 


1.46814 


0.47750 


2.53275 


1.45136 


0.86548 



Table 1: Numerical data for p = 0, 0.1 , 0.2 , 0.3, 0.4. 



Another important problem is to investigate the relation between the flux 
and the density, which is called fundamental diagram. The density p of the 
cars is defined by N/L, where we choose L=200 and vary from 10 and 300 
in the simulation. The flux Q is defined by the number of cars passing by a 
position per unit time. The data was accumulated and averaged over during 
20,000 time after first 1000 time. Numerical results are plotted in Figure | 
for p = 0, 0.1, 0.2, 0.3. Here we omitted to plot the data which overwrites 
the other data. 
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Figure 4: Flux-density fundamental diagram for p = 0, 0.1 , 0.2 , 0.3. 
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In the homogeneous flow, the relation between flux Q and density p is 
given by 

Q = pV = p (^tanh Q - 2^ + tanh(2)^ . (17) 

In Figure it is represented by the dashed curved line. The numerical results 
agree with this line in the homogeneous flow region. 

As is discussed in Ref. 0, the flux-density relation in the congested flow is 
evaluated from the data in Table |l|. In the congested flow region, the relation 
between the density p and the flux Q is given by 



Q = - HackP. (18) 

Substituting the values in Table |1| into eg. (p!8|) , we obtain the Q-p relation 
as in the Table ^. These lines are plotted in Figure ^ with the fundamental 



p 


Q-p relation 


0.0 


g = 0.55597- 0.14792 p 


0.1 


Q = 0.63853 - 0.31302 p 


0.2 


Q = 0.73174 - 0.49945 p 


0.3 


Q = 0.82518 - 0.68632 p 


0.4 


Q = 0.91475 - 0.86548 p 



Table 2: Flux-density relation in congested flow for p = 0, 0.1, 0.2, 0.3, 0.4. 



diagram. In Figure 4 the vertical lines represent the boundaries for the 
stability condition given by eq.(|l^) in the case of p = 0, 0.1, 0.2, 0.3, where 
the same dashed line is used as the one used to draw the predicted line in 
eg . (plSf) . The numerical results for various p values are good agreement with 
the predicted lines in Table 0. 

We will summarize the effect of the Ax„+i dependent term in fundamental 
diagram. In the homogeneous flow region, there is no effect in the flux-density 
relation, because p dependence in the generalized OV function disappears in 
the case of homogeneous flow, as is seen in eq.(^. In the congested region, 
the flow increases as p increases if p < | and inversely the flow decrease as 
p increase if p > |. The reason is as follows. When the density is low, the 
average of the headway is long. If Axn+i is long, the larger value of the 
velocity than that of the case without Axn+i dependence is allowed. Hence 
the flow increases by taking into account the headway of the immediately 
preceding car. Inversely, when the density is high, i.e. p > ^, the small value 
of the velocity is taken compared with the case of p = 0, because Axn+i is 
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short. Hence the flow decreases by taking into account the headway of the 
immediately preceding car. 

4 Rescaled Model 

As far as the stability condition of traffic flow is concerned, our analytical 
result of eq.(|13[) means that the effect of the Axn+i dependent term can 
be compensated by rescaling the sensitivity a. Thus we can rewrite the 
generalized OV model as 

Now the stability condition of the above model is given by / < | which is 
independent of p. By investigating the hysteresis loop in this model with 
various values of p, we can clarify the effect of the headway of the preceding 
car. Numerical simulation can be performed in the same way as in Section 
3. Obtained numerical data for the headway and velocity of a particular car 
is plotted in Figure |. 




0.5 1 1.5 2 2.5 3 3.5 4 
headway A 

Figure 5: Hysteresis loops for p = and p = 0.2. 

We can find that the effect of the Axn+i dependent term can not be 
compensated by rescaling the sensitivity a. The effect of p is now easily 
understood. Compared with p = 0, shorter value of the minimum of the 
headway is allowed, and also longer value of the maximum of the headway is 
allowed. Furthermore, when the car accelerates, larger value of the velocity 
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can be taken even if the headway is short. Inversely, when the car decelerates, 
smaller optimal velocity is taken even if the headway is long enough. 

The above characteristic feature in the phase space holds for various p 
values as long as p < |. If we take p > ^, the numerical simulation shows 
that the cars take over the cars ahead, thus it is not realistic. Of course, a 
room to change the form of OV function is left, but we will not consider it 
here. 

It is found that the model give by eg. (pJ]) has a different feature in the 
unstable state when we take p = ^- Figure || shows the hysteresis loop for 
P = ^ after some relaxation time, e.g. about 10*^ time in the simulation. 
More precisely, two independent hysteresis loops appear for two successive 
cars. In Figure ^, the line and the dotted line are the hysteresis loops for the 
odd car number and even car number, respectively. This hysteresis loop is 
well understood when the orbit of the car is shown in the (Ax^, Ax„+i, x^) 
phase space. Figure shows the orbit of two successive cars in the phase 




0.5 1 1.5 2 2.5 3 3.5 4 
Headway A X^, 



Figure 6: Hysteresis loops for p = ^ Figure 7: Hysteresis loops for p = ^ 
in (Ax„,i„) space. in (Ax^, Ax„+i, Xn) space. 

space. The simulation shows that all of the vehicles sweep the same orbit 
as the one in Figures |^ and |^. We have examined the model with various 
initial configurations such as random distributions, different number of cars 

and circuit length L. The above result does not depend on the initial 
distributions and the parameters. 

Figure | and ^ show the headways and the velocities for all vehicles, re- 
spectively, in the case of = 100 and L = 200 at 10^ time, starting with 
arbitrary random distributions with zero velocities. Numerical simulations 
show that the characteristic feature does not depend on the initial distribu- 
tions and the parameters. We can find that an intermediate state appears 
for the movement of cars, which is different from the ordinary stop-and-go 
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Figure 8: Headways for all vehicles Figure 9: Velocities for all vehicles 
at 10^ time with N = 100, L = 200. at 10^ time with N = 100, L = 200. 



state. Numerically, the state corresponds to the points (Ax„, Ax„+i, x„) = 
(3.677, 0.323, 0.964) and (Ax„, Ax„+i, x„) = (0.323, 3.677, 0.964) in the 
phase space as is given in Figure |^. Thus the intermediate state means 
that the car can keep a certain velocity (0.964) though the headway is short 
(0.323) or long (3.677) with holding Ax„ + Ax„+i = 4. The appearance of 
the intermediate state is universal in the sense that it does not depend on 
the initial distributions, the number of cars and the length L. It should 
be noted that this state never appears if p < ^. 

This characteristic feature is understood by considering the equation for 
the headway which is given by 

= - 2p)V{AXn+l) + pV{AXn+2) - (1 - p)V{AXn) - AXn) ■ 

l + 2p 

(20) 

If p = I is chosen, the first term in the right hand side of eq.(20) vanishes and 



hence Axn+i dependence disappears. One might expect that the intermediate 
state is understood by the nonlinear analysis near the critical point a = ac- 
Following the analysis by Komatsu and Sasa[^, the equation for the headway 
is derived near the critical point and the modified Korteweg-de Vries(MKdV) 
equation and its higher-oder corrections is obtained by introducing a small 
scaling parameter e = ^y {ac — a) /etc- However, the equation near the critical 
point in the case of p = | is exactly the same as the one in the case of p = 
up to the order e. Because the difference in eq. (|20D with p = and P = | 
near the critical point is that the Fourier mode with p = | is just twice that 
of the other one with p = 0. Since our analysis is performed far from the 
critical point, the higher-oder corrections or non-perturbative effect should 
be considered to understand the intermediate state. 
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5 Summary and Discussion 



We have analyzed the stabihty of the generahzed optimal velocity model 
where the optimal velocity function depends not only on the headway of each 
car but also on the headway of the immediately preceding one. The effect 
of the newly introduced Ax„+i dependent term was examined by numerical 
simulation. In particular, the hysteresis loop in the phase space and the 
flux-density relation were examined in detail by taking the various values of 
the parameter. 

We found that the effect of the Ax„+i dependent term can not be com- 
pensated by rescaling the sensitivity a. In the model with the specific pa- 
rameter choice p = |, we found that the intermediate state appears for the 
movement of cars, which is different from the ordinary stop-and-go state. 
Numerical simulation shows that the appearance of the intermediate state is 
universal because it does not depend on the initial conditions. 

We would expect that the model is related with the exact solution given 
by Jacobi's elliptic function [|TO|] . It is interesting to examine the difference- 
differential equation 

Xn{t + t) = V{AXn{t), Ax„+i(t)). (21) 

The details of the numerical simulation and the analysis of the underlining 
mathematical structure of the generalized optimal velocity model are under 
study. 
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